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Abstract 

This is an elementary introduction to the recent significant theoretical progress in the field 
of heavy quarkonium physics. We show how renormalon cancellation takes place in the 
heavy quarkonium system, such as bottomonium and (remnant of) toponium resonance, 
and how this notion is useful in extracting the MS masses of the bottom and top quarks. 



* Based on the invited talk "Top quark physics at future linear colliders" given at the Japan Physics Society 
Meeting, Osaka, Japan, March 30 - April 2, 2000. 



1 Introduction 



Recently there has been significant progress in our understanding of heavy quarkonia such 
as T's and remnant of toponium resonances. Developments in technologies of higher order 
calculations and the subsequent discovery of renormalon cancellation enabled extractions of 
rrib and (in future experiments) of rrit with high accuracy from the quarkonium spectra. In 
this article we review the notion of renormalon and its cancellation in the heavy quarkonium 
system. We demonstrate how it is useful in extracting the quark masses.* 

We consider heavy quarkonia whose sizes (given by the Bohr radius ~ (a^w^g)^^) are much 
smaller than the hadronization scale Aqqj-, ~ (0.3 GeV)~^. In reality the candidates are T(15') 
and (remnant of) toponium resonances, whose sizes are ~ (1.5 GeV)~^ and ~ (20 GeV)~^, 
respectively. In such a system, gluons participating in the binding of the qq boundstate have 
wavelengths much shorter than the hadronization scale, so theoretically nature of the bound- 
state can be described well using perturbative QCD. In particular the boundstate spectrum (the 
mass of boundstate) can be calculated as a function of the quark mass and as- Consequently 
we can extract the quark masses, m^, m^, from the masses of the above quarkonia. 

First let us state briefly the theoretical framework used in contemporary calculations of 
spectra of non-relativistic boundstates such as the heavy quarkonia. In old days people solved 
the celebrated Bethe-Salpeter equation to compute the boundstate spectrum. We no longer 
use this equation; instead we reduce the problem to a quantum mechanical one. Namely we 
solve the non-relativistic Schrodinger equation 



to determine the boundstate wave functions and energy spectrum. The quantum mechanical 
Hamiltonian is determined from perturbative QCD order by order in expansion in 1/c (inverse 
of the speed of light): 



Since quark and antiquark inside the heavy quarkonium are non-relativistic, the expansion in 
1/c leads to a reasonable systematic approximation. Presently the Hamiltonian is known up to 
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* See Ref. for a comprehensive review of renormalo: 
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where m denotes the pole mass of the quark; as = a;s(yu); Cp = 4/3, Ca = 3 are color factors; 
fi' = fi e'^^ . The lowest-order Hamiltonian Hq is nothing but that of two equal- mass particles 
interacting via the Coulomb potential. 

In addition to the above Hamiltonians, some of the terms in higher order Hamiltonians if„ 
(corresponding to the underlined terms) are known. From the analysis of these higher order 
terms, one finds that there exists a problem in extracting the quark mass from the boundstate 
spectrum. We will first address the problem, which is known as the "renormalon problem", 
and then will see how it is solved. 



2 The Renormalon Problem 

The terms with underlines in Eqs. stem from the running of the coupling constant 

dictated by 

and all higher order terms can be determined using the renormalization-group equation. Thus, 
in the "large (3q approximation" , the potential between quark and antiquark is given by[| 

^/3oW = -y ^e'? C'i. ; g=|gl, (7) 

where the 1-loop running coupling is defined as a perturbation series in ^^(/i): 

a.(,)£{-^.og(4)}". (8) 

n=0 \/i / J 

We may not resum the geometrical series before the Fourier integration because then the inte- 
grand exhibits a pole at g = A, 

^ - ^'/^° ^ (9) 



and the Fourier integral becomes ill-defined. Here, 

A = /i exp 



(10) 



is the /i-independent integration constant of the 1-loop renormalization-group equation. There- 
fore, the above potential Vp^ir) can only be defined as a perturbation series in 0:5 (/x). (Fourier 
integral of each term of the series is well-defined.) 

T It is the Coulomb potential with the "running charge"; of. Eq. (|). 
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When we examine the large-order behavior of this perturbation series, we find that it is an 
asymptotic series and has an intrinsic uncertainty 



(5l^^,(r) ~ A = C(300 MeV). (11) 

If we want to extract the quark mass from the spectrum of boundstates, this uncertainty in the 
potential is directly refiected to the uncertainty in the quark mass. It is because the quarkonium 
mass is determined as twice the quark pole mass minus the binding energy, and an uncertainty 
in the potential means an uncertainty in the binding energy. This implies that we cannot 
determine the quark mass to an accuracy better than (9(300 MeV). 

Now we examine the series expansion of Vpglr) and see how the uncertainty arises ^. If 
we perform the Fourier integration term by term, 

= -CFAnasifi)f:\^^^^y fn{r,ix)xn\, (13) 
n=0 ^ 47r J 

the coefficients fn{r,fi) can be determined from a generating function 

nr,,,n) . f^^(^ (14) 



(27r)3 g2 ^^2 
(/ir)2" 



(15) 



Anr cos(7ru) r(l + 2u) 
= ^/„(r,/x)«". (16) 

n 

Using this generating function, one easily obtains the asymptotic behavior of fn{r,fi) for large 
n. The large-n behavior of /„(?", A'-) determines the domain of convergence of the series expansion 
(p!6|) at M = 0. Conversely from the structure of the pole of ([T5|) nearest to m = (see Fig. |l|), 
which limits the radius of convergence, one obtains^ for n ^ 1 

/„(r,/i) ~ ^/^x2". (17) 

Note that this asymptotic behavior is independent of r. This means that, although each term of 
the potential Va(,(r) is a function of r, its dominant part for n ^ 1 is only a constant potential 
which mimics the role of the quark mass in the determination of the quarkonium spectrum. 
Thus, asymptotically, the n-th term of V^^, (r) is given by 



The leading asymptotic behavior of /„(r, /i) is same as that of the expansion coefficients of Res[-F'; u 
i]x(«-i)-i. 
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Figure 1: Analyticity of the generating function F{r,fi;u) shown on the complex w-plane. Poles are located 



at u = 



13 5 
2 ' 2 ' 2 ' 



•. Also the domain of convergence of the series expansion at u = is shown. 




Figure 2: The graph showing schematically the asymptotic behavior of the n-th term of —Vf3g{r) for n 2> 1. 



As we raise n, first the term decreases due to powers of tfie small as', for large n the term 
increases due to the factorial n\. Around no = 27r/{l3oas{n)), V^^^ becomes smallest. The size 
of the term scarcely changes within the range n e (tzq — yno, uq + ^/no); see Fig. ^ We may 
consider the uncertainty of this asymptotic series as the sum of the terms within this range: 



6Vf 



r] ~ 



E 

n=no-^/nQ 



V,. 



(n) 



A. 



(19) 



The //-dependence vanishes in this sum, and this leads to the claimed uncertainty. 

In passing, we note that this asymptotic series is not Borel summable; a Borel summable 
series has terms alternating in sign, but the asymptotic series originating from the QCD infrared 
renormalon has terms with the same sign. We cannot circumvent the uncertainty by Borel 
summation of the series. 
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Figure 3: The 1-loop quark self-energy diagram and the one-gluon-exchange diagram between quark and anti- 
quark, where the gluon propagators are replaced by the bubble- chains. These diagrams, respectively, contribute 
to the renormalons in 2mpoic and V/}„{r). 



3 Renormalon Cancellation in the Total Energy of a qq 
system 

Now we state how the problem can be circumvented. Consider the total energy of a color-singlet 
non-relativistic quark-antiquark pair: 



It was found that the leading renormalon contained in the potential V^o(r) gets cancelled 
in the total energy Etot{f) if the pole mass mpoie is expressed in terms of the MS mass. The 
potential and the pole mass are expressed in terms of the 1-loop running coupling aii^{q) as 



The potential Vf3g{r) is essentially the Fourier transform of the Coulomb gluon propagator 
exchanged between quark and antiquark; the difference of mpoie and mjj^ is essentially the 
infrared portion of the quark self-energy, see Fig. ^. As we saw, the renormalon uncertainty is 
related to the "would-be pole" contained in aii^{q), cf. Eq. The signs of the renormalon 
contributions are opposite between Vgo(^) ^poie because the color charges are opposite 
between quark and antiquark while the self-enregy is proportional to the square of a same 
charge. Their magnitudes differ by a factor of two because both the quark and antiquark 
propagator poles contribute in the calculation of the potential whereas only one of the two 
contributes in the calculation of the self-energy. Expanding the Fourier factor in 2mpoie in a 
Taylor series for small q, 



\otir) ^ 2mpoie + Vao(r). 



(20) 




(21) 



(22) 



e 



l+iq- f+ -{iq - rf' + ---, 



(23) 
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the would-be pole contained in the leading term gets cancelled against Vao(^)| consequently 
the renormalon contributions cancel. 

As a result of this cancellation, the series expansion of the total energy in asifJ') behaves 
better if we use the MS mass instead of the pole mass. Residual uncertainty due to uncancelled 
pole can be estimated similarly as in the previous section and is suppressed as 




which is much smaller than the original uncertainty. 



4 Extracting the MS masses using the Full NNLO Result 

To see how well the renormalon cancellation works, we examine extractions of the bottom and 
top quark masses using the full next-to-next-to-leading order (NNLO) result of the boundstate 
spectrum. The full NNLO formula for the lowest lying (IS*) boundstate can be calculated from 
the Hamiltonian Eqs. (^-(^ and is given by 0, 



4 

Mis = 2mpoie - -asififnipoic 
y 



I ( 1793 I 29177r^ _ ^ t 275C3 A , / 1693 _ Utt^ _ 19C3A„ t { JL eL t 2(3 A ^2 
"•"V 12 216 32 4 y ^ V 72 18 2 y""' V 108 54 9 J "'I 



(25) 



where L = log[/i/(Ciras'(/i)mpoic)], and rii denotes the number of massless quarks. We may 
examine the size of each term of the above perturbation series. Alternatively we may rewrite 
the above expression in terms of the MS mass and examine the series. Presently the relation 
between mpoie and m = '"^ms("^ms) known up to three-loop order p: 



mpoie = m X 



1 + -(^^) + [^^)\-lMUni + 13.4434) 



+ (^^^)^(0.6527nf - 26.655n« + 190.595) 



(26) 



First we apply the formula to the T(15') state, for which rii = 4. Taking the input parameter 
as mpoie = 4.96 GeV/m = 4.22 GeV and setting fi = m (i.e. expansion parameter is asijn) = 
0.22), 

Mt(is) = 2 X (4.96 - 0.05 - 0.08 - 0.11) GeV (Pole-mass scheme) (27) 
= 2 X (4.22 + 0.35 + 0.12 + 0.04) GeV (MS-scheme). (28) 



^ We are interested in the infrared region q ^ A. The expansion is justified since the typcal distance between 
quark and antiquark is much smaller than the hadronization scale (Ar) <^ 1. 
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Figure 4: The total production cross sections for e+e^ tt in the threshold region, where the leading-order 
(LO), next-to-leading-order (NLO) and next-to-next-to-leading-order (NNLO) curves are shown. The vertical 
lines represent the locations of the corresponding IS" resonances when the top width is artificially taken to zero; 
also the position of 2mt is shown by vertical lines. The two figures correspond to (a) the pole-mass scheme, and 
(b) the MS-scheme. These figures are made by T. Nagano. 



One sees that the series is not at all converging in the pole-mass scheme, whereas in the MS- 
scheme the series is converging quite nicely up to the calculated order. (See Sec. || for details 
of how we derived the series in the MS-scheme.) Comparing this with the experimental value 
^T(is) = 9.46037 ± 0.00021 GeV, one may extract the MS bottom quark mass 

nib = "^Ms(^Ms) = 4-22 ± 0.08 GeV. (29) 

One might think that, looking at the behavior of the above series, we may assign a smaller 
theoretical uncertainty. The present uncertainty is, however, dominated by non-perturbative 
uncertainties other than the renormalon contributions. Thus, presently is determined to 2% 
accuracy [§, |T^. It seems to be fairly good in view of the fact that its major part is controlled 
by perturbative QCD. 

Next we turn to the (remnant of) "toponium", for which rii = 5. At future linear e~^e~ or 
colliders, the top quark mass will be determined to high accuracy from the shape of the 
tt total production cross section in the threshold region. The location of a sharp rise of the 
cross section is determined mainly from the mass of the lowest lying (IS*) tt resonance, so we 
will be able to measure the resonance mass and extract the top quark mass. Similarly to the 
previous case, we set mpoie = 174.79 GeV /m = 165.00 GeV, asijn) = 0.1091 and obtain 

Mis = 2 X (174.79 - 0.46 - 0.40 - 0.28) GeV (Pole-mass scheme) (30) 
= 2 X (165.00 + 7.20 + 1.24 + 0.22) GeV (MS-scheme). (31) 

In Figs. § are shown the convergence properties of the above series together with the correspond- 
ing cross sections. In the pole-mass scheme, the convergence is very slow. According to the 
renormalon argument, also uncalculated higher order terms would not become much smaller. 
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On the other hand, in the MS-scheme the series shows a healthy convergence behavior. For top 
quark, non-perturbative uncertainties are much smaller than the present perturbative theoreti- 
cal uncertainty. Thus, from the above series we estimate that can be determined to around 



100 MeV accuracy [11 



5 Physical Implications 

Let us discuss some physical implications of renormalon cancellation in the heavy quarkonium 
system. Firstly, as already mentioned, we expect that gluons with wavelength much longer 
than the size of the quarkonium cannot couple to this system. Hence, we expect that infrared 
gluons with momentum transfer q <C asm should decouple from the expression of -Etot- This 
is a naive expectation based on classical dynamics. Such understanding should be valid when 
it is described by the bare QCD Lagrangian without large quantum corrections. The MS mass 
is closely related to the bare mass of quark; only ultraviolet divergences are subtracted. On 
the other hand, the pole mass has much more intricate relation to the bare mass, because the 
relation includes in addition infrared dynamics of the quantum correction to the quark self- 
energy. In this sense it would be natural to expect the decoupling phenomenon to be realized 
when the MS mass is used to express i^tot- 

Secondly, the pole mass of a quark is ill-defined beyond perturbation theory. It can be 
determined only when the quark can propagate an infinite distance. Generally accepted belief 
is that when quark and antiquark are separated beyond a distance ~ the color flux is 
spanned between the two charges due to non-perturbative effects and the free quark picture is 
no longer valid. On the other hand, it is natural to consider the total energy (or the mass) of 
a quarkonium which is a color-singlet state. It can propagate for a long time and the notion of 
its mass is not limited by the hadronization scale. 

Thirdly, renormalon cancellation seems to be a universal feature which occurs process inde- 
pendently. The same phenomenon was known e.g. in the QCD corrections to the p-parameter 
||r^ and B decays [T^. Therefore, the MS mass, which is determined accurately from the 
quarkonium spectum, would be more suited than the pole mass for an input parameter in 
describing other physical processes. 



6 How to Cancel Renormalons in the Quarkonium Spec- 
trum 

There is one non-trivial point in realizing renormalon cancellation in the perturbation series 
of the quarkonium spectrum. When the pole mass and the binding energy are given as series 
in as-, renormalon cancellation takes place between the terms whose orders in as differ by one 
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Figure 5: The figure showing how cancellations should take place between the diagrams. The orders of the 
potential graphs are shifted by one power of as which is provided by the inverse of the Bohr radius ( 7) ~ asm. 
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2m 




B2 al + 53 a% 



(32) 



cf. Eqs. (P3| ) and (^). Intuitively this can be seen from the diagrams shown in Fig. |^. An 
additional power of as in the binding energy is provided by the inverse of the Bohr radius 
(r~^) ~ asm, 



I.e. 



ai 



(33) 



Still, one might wonder how cancellation can ever take place between different orders in as 
for any value of as- So we demonstrate the cancellations at large orders in a specific example. 
The mass of the IS* boundstate can be written in the form of an expectation value as 



M,s = {lS\2m^,,, + H\lS) 



(34) 



using the Hamiltonian (||) and its IS energy eigenstate 1 15"). We are interested in the leading 
renormalon contributions, so we replace 



2mpoie + H 



2m 



pole 



(35) 
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The energy eigenstate can be expanded in 1/c: 

\is) = \ is^^^) + -\is^'^) + his^'^) 



(36) 



Renormalon cancellation takes place in arbitrary combination of (IS'*'*-'! ■ ■ ■ llS'*--')), but for sim- 
plicity we evaluate only the following part: 



(0) 
IS 



(37) 



The second term corresponds to the binding energy and we may evaluate it at each order of 
the perturbation series: 



(15W|r,„(r)|15(°)) 



drr^ |i?(r)p\/^,(r) 
-CpasW mpoic X 2^ <j — — — I gn{nao)xnl, 



n=0 



where the zeroth-order IS Coulomb wave function is given by 



R{r) 



3/2 



-r/ao 



ao 



^Cpas{fi)mpo\c ) : Bohr radius. 



(3^ 



(39) 



5'„(/iao)'s are polynomials of log(/iao). Using the generating function method, one obtains the 
asymptotic form 

2 



5(„(/iao) ~ -/iao X 2" oc 

IT as 



(40) 



Thus, for n ^ 1, it becomes proportional to a^^ and effectively shifts the order of as- By 
setting fi = m, it is easy to check that in this example the leading renormalon cancels as in 
Eq. ( ^2]) and the residual piece behaves as 

' Poasim) 



n-th term of M-, 



(0) 
IS 



asym) m x n\ x 



It follows that 



A X 



\asm/ 



< A. 



(41) 



(42) 



From this example one learns that the cancellation of renormalon contribution between 
shifted orders An as and -B„+i a^"*"^ should be properly taken into account when expressing the 
boundstate mass as a perturbation series. There are many different prescriptions to accomplish 
this. We derived the series (|2^) and (^T]) in the following manner. We have rewritten Eq. (^) 

as 



Ml 



IS 



2m 



pole 



2m X 



X 



1 + PnOis{m 

n=2 

3 

1 + ^ Anasim) 
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n=l 



1 + ^ Pnas[ 



m] 



n=2 



(43) 



where P„'s are polynomials of \og[as{m)j and A^s are just constants independent of 0:5(771). 

We identified PnO^s order a^"^ and then reduced the last line to a single series in 05. 

Parametric accuracy of the last terms of Eqs. (^) and (|3T|) is a|m. In order to improve 
the accuracy to agfn, we need to know further (a) the exact 4-loop relation between ?7ipoie and 
m, and (b) the binding energy at order a| in the large-/3o approximation. 

7 Some Questions 

One can ask some interesting questions related to the renormalon problem and extraction of 
quark mass. In the case of top quark, its mass will also be measured from the invariant mass 
distribution of decay products of the top quark in future hadron collider experiments and in 
future e~^e~ collider experiments. What is the mass extracted from the peak position of the 
Breit-Wigner distribution? Naively it is the pole mass because the peak position is determined 
from the position of the pole of the top quark propagator in perturbative QCD. On the other 
hand, as we have seen, the pole mass suffers from a theoretical uncertainty of (9(300 MeV). 
Since future experiments will be able to determine the peak position of the invariant mass 
distribution to 0(100 MeV) accuracy, we will indeed face this serious conceptual problem. 
We believe that renormalon cancellation also takes place in this physical quantity. But the 
problem lies in the fact that we have no reliable theoretical method to calculate the invariant 
mass distribution of realistic color-singlet final states. Rather we calculate the invariant mass 
of final partons which do not combine to color-singlet state. Another question is whether the 
peak position of invariant mass distribution measured in hadron collider experiments will be 
the same as that measured in e^e~ collider experiments. The answer is probably no, since at 
hadron colliders top and antitop are pair-created not necessarily in a color-singlet state. 
These are very interesting questions which are worth further studies. 
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